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Abstract 


A theory is presented for computing waves radiated from waves 
on a smooth surface. With the assumption that attenuation of the 
surface wave is due only to radiation and not to dissipation in the 
surface material, the radiation coefficient is derived in terms of the 
attenuation factor. The excitation coefficient is determined by the 
reciprocity condition. Form ulas for the shape and the spreading of the 
radiated wave are derived, and some sample calc illations are pre sented. 
An investigation of resonant phase matching for nonseparable surfaces 
is presented with a sample calculation. A discussion of how such 
calc illations might be related to resonant frequencies of nonseparable 
thin shell structures is included. A description is given of nonseparable 
surfaces that can be m odeled in the vector form that facilitates use of 
the appropriate form ulas of differential geom etry. 


1. Introduction 

Calculations of surface waves and of waves radi- 
ated from surfaces have potential applications in such 
areas as sound diffraction and shock diffraction by an 
aerodynamic surface, nondestructive testing (crack 
detection), mutual coupling of surface-mounted an- 
tennas, underwater signal recognition, and resonance 
of thin shell structures. Thg geometric theory of 
diffraction (GTD) developed by J. B. Keller and his 
associates (ref. 1) three decades ago led to many 
such applications and to theories for diffraction by 
a smooth body (ref. 2) and for surface waves on a 
curved surface (ref. 3). This development, applica- 
ble as a high frequency approximation, represented a 
significant breakthrough in the wave theory involving 
surfaces for which the wave equation is not separa- 
ble. By the GTD method, the surface waves could be 
calculated by ray tracing techniques; then related to 
waves in the surrounding medium by means of coef- 
ficients, variously termed diffraction, excitation, ra- 
diation, or launching coefficients. These coefficients 
could be determined by exact solutions to canonical 
problems, that is, tractable problems involving sepa- 
rable surfaces. Thus, a large class of smooth surfaces 
became amenable to surface wave analysis. 

However, the results, in terms of actual use, have 
been disappointing. A cursory survey of the litera- 
ture yields few, if any, actual GTD calculations in- 
volving nonseparable surfaces. Reference 4 describes 
in detail a procedure for calculating ray paths and 
wave fronts on such a surface. The purpose of this 
paper is to supplement reference 4 by relating the sur- 
face waves to (1) the diffracted waves or head waves 
radiated into the surrounding medium, and (2) nat- 
ural resonance phenomena of thin shell structures. 
Equations for the radiated waves are derived, and 


sample calculations illustrate the procedure. (Head 
waves are waves radiated at an angle not equal to 
0 °.) 

Rays that traverse the surface along a closed path 
give rise to a resonance effect at critical frequencies 
(ref. 5). This effect is observed in the diffracted or 
head wave radiated from the surface. A method for 
locating closed paths associated with this kind of 
resonance is described, and an illustrative calculation 
is included. 

2. Symbols 

A amplitude 

AA wave front vector, pointing in 

direction of wave front normal, with 
magnitude equal to ray tube cross- 
sectional area 

c m wave speed in medium surrounding 

surface 

c.< 5 , surface wave speed 

E, F, G metric coefficients 

E r , G r metric coefficients along and normal 

to surface ray 

A E energy flux 

e,f,g coefficients of second fundamental 

form (curvature coefficients) 

I incident ray direction 

I|| vector projection of I onto tangent 

pi ane 

i,j,k orthonormal base vector in direction 

of x , y, and z axis, respectively 
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total (Gaussian) curvature 
wave number 

number of times surface ray tra- 
verses a caustic 

unit surface normal vector 

surface distance normal to ray 
direction 

unit vector defined by equation (32) 

unit vector normal to p in T , N 
plane 

radiation coefficient 

vector location of generic surface 
point 

distance along ray in surrounding 
medium 

distance along ray from surface to 
point P 

distance along ray from point P to 
surface 

unit tangent vector to surface ray 
time on surface ray 
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Subscripts : 
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P, Q, P r , Q j 

r 

s 


curvature of surface ray 

parameter defining a local direction 
on surface (see eq. (17)) 

attenuation factor 
homotopy function 
scaling, or sizing, function 
density 

distance along ray path 
cross-sectional area of ray tube 
torsion of surface ray 
ph ase 

frequ enc y, rad ians 

initial 
fin al 

condition in surrounding medium 
initial point of surface ray 
quantity evaluated at this location 
surface ray 
surface 


time on ray emitted from surface 


3. General Considerations 


surface coordinates 

unit vector normal to p and q 

Cartesian coordinates 

vector location of point on surface 
ray 

vector location of point on ray 
emitted from surface 

coordinates of points defining 
normalized cross-section shape 

excitation coefficient 

parameter distinguishing individual 
rays emitted from surface at time t 

angle of incidence or radiation 


First, section 4.1 describes the basic assumptions 
and scope of the theory. Then section 4.2 demon- 
strates how a general class of smooth surfaces can be 
modeled in an appropriate analytic form for the ap- 
plication of the theory and procedures described in 
the subsequent sections. This section is a brief sum- 
mary of the method of reference 6. Section 4.3 dis- 
cusses calculation of surface waves and how they are 
related to waves in the surrounding medium. Sec- 
tion 4.4 is a discussion of related problems. The 
strengths of the radiated waves as well as the global 
representation of their shapes are treated in sec- 
tion 4.5. Several sample calculations of these wave 
shapes are included in section 5. Section 6 discusses 
the theory and calculation procedure for finding i In; 
closed surface ray paths that give rise to resonance 
effects. 


geodesic curvature of surface wave 
front 

normal curvature of surface wave 
front in direction p 

curvature of surface wave front in 
direction p 


4. Analysis 

4.1. Scope of Theory 

Throughout the analysis, the emphasis is placed 
on theory and methods for determining the geometric 
quantities that are required for computing wave fields 
that are emitted by waves on a smooth surface. The 
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primary assumptions and limitations of the theory 
are (1) scalar waves, (2) real rays and surfaces, 
(3) homogeneous surfaces, (4) surrounding medium 
can support a wave, and (5) radiated rays do not 
intersect surface rays. An attempt has been made 
to minimize, the complications associated with the 
specialized physics associated with different types of 
waves. For example, the discussion is confined to 
scalar waves. However, in much of the literature 
on vector waves (ref. 7, for example) the results are 
synthesized from combinations of results for scalar 
waves. A brief discussion is included in section 4.4 on 
some additions and modifications that are required 
for some particular applications. 

The elegant concept of complex rays, utilized in 
reference 3, is not used in this analysis. Although 
mathematically efficient, such concepts tend to in- 
timidate many engineers. Only real rays, real sur- 
faces, and real frequencies are treated, with damping 
due to radiation considered as attenuation and not 
as an imaginary part of the phase. 

Derivations of surface wave speeds have been 
given in a multitude of papers (ref. 8 and its ref- 
erences) and is not duplicated here. It is assumed 
that the surface wave speed is known and that it is 
constant (the surface is homogeneous). It should not 
be too difficult to extend the analysis and methods 
to inhomogeneous surfaces inasmuch as a procedure 
for calculating surface ray paths on such surfaces has 
already been described (ref. 4). However, the anal- 
ysis is considerably shortened and simplified by re- 
stricting it to homogeneous surfaces for which the 
ray paths are geodesics. 

The surface is assumed to be immersed in a 
medium that will support a wave whose speed is 
equal to or less than that of the surface wave speed. 

Finally, assumptions concerning the convexity of 
the surface need to be considered. The usual re- 
striction is to smooth convex surfaces. Requiring 
smoothness avoids the problems associated with the 
specialized treatments required for diffraction ef- 
fects resulting from discontinuities in slope or curva- 
ture,; Requiring convexity avoids the complications 
that arise when the surface waves exist simultane- 
ously with incident waves, either from the off-surface 
source (fig. 1(a)) or from previously diffracted waves 
from the surface itself if it is inflected (fig. 1(b)). 
However, a global requirement of convexity is often 
too severe h restriction. The only restriction required 
in order to secure the desired simplification is that, 
for any ray under consideration, the normal curva- 
ture must be positive outward over that part of the 
ray path that is being investigated. For example, if 


the object is illuminated from a far-held source so 
that the incident rays all have essentially the same 
direction, and if the grazing rays diffract along paths 
that satisfy this requirement, then no difficulties oc- 
cur if the surface is concave in the direction Orthogo- 
nal to the surface rays (top part of fig. 2). As a second 
example, the problem of mutual coupling of surface- 
mounted antennas depends only on the geodesic ray 
path connecting the two antennas. Consequently, 
this problem can be treated if the ray path lies along 
a convex part of the surface. As a final example, if 
one is considering the head waves excited by surface 
waves which move at a speed several times gf later 
than the wave speed in the surrounding medium (as 
is usually the case) then the requirement for convex- 
ity, even in the ray direction, can be relaxed except 
in the most severe cases. For these waves, the rays 
are radiated from the surface at such a large angle 
that they are not likely to strike another part, of the 
surface (bottom part, of fig. 2). 



(a.) Interaction with incident rays. 



(b) Interaction with diffracted rays. 

Figure 1. Interaction of surface rays with incident or diffrac- 
ted rays on inflected surface. 

4.2. Analytically Lofted Surfaces 

As was stated in section 1, the analysis and pro- 
cedures to be described are applicable to a. general 
class of smooth surfaces and are not. restricted to 
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Figure 2. Problems for which theory is applicable on inflected 
surfaces. 


surfaces for which the wave equation is separable. 
The analysis assumes only that the surface can be 
expressed as a vector function r of two variables: 


the cross sections can be specified by prescribing a 
scaling function v(x), which varies smoothly but not 
necessarily monotonically from an initial to a final 
value as x varies from x = x; to x = x f. The scaled 
variables are defined by 

y(x,£) = v{x) y(x,£) (3a) 

r(a;,£) = v{x) z(x ,£) (3b) 

The vector equation of this transition surface is 

r = xi + //I -c . £ !,/ + z(x,£)k (4) 

which has the required form (eq. (1)) where the sur- 
face variables (u,v) are (x,£). For modeling com- 
plicated shapes, v may be allowed to vary with £ as 
well as with x, or one may take different scaling func- 
tions I I'y. /'■ ! for the (y,z) coordinates. However, for 
illustrative purposes, equations (3) describe a suffi- 
ciently general class of surfaces. The x = Constant 
lines describe the cross-section shapes,,, and the £ = 
Constant shapes describe the lofting lines. The first 
derivatives are 

\dv . d u , 

,v = ,+ + J 


r = r ( u ,v) ( 1 ) 

A class of such surfaces, introduced in reference 6 
for the purpose of expediting design of aerodynamic 
components, is described here. These surfaces are 
defined by specifying the shape of an initial cross 
section in a plane x = x, and a final, or terminating, 
shape at x = x f and then requiring the intermediate 
cross sections to represent a smooth transition from 
the initial to the final shape. The shapes (t/(£), -(£)) 
of the two extreme cross sections are described para- 
metrically in terms of an arbitrary parameter. This 
parameter might be, for example, the normalized y 
coordinate, the normalized arc length, the angle in 
polar coordinates, or the ellipse angle parameter. 

Let y{x) be a function, varying smoothly and 
monotonically from 0 to 1 as x varies from x, to xf 
(a homotop y function). Then, for Xj < & < xj-, a 
cross-section shape is defined by 

y(x, £) = [1 - y(x)\ */;(£) + y(x) y f (£) (2a) 

i( -c, £) = [1 - y(x)] + y(x) zf(£) (2b) 

These shapes vary smoothly from the initial shape at 
x = x - to the terminating shape at x = x j - . By taking 
£ as a normalized coordinate, the size variation of 


I" dv dy 1 

+ (5a) 

M n d Vi _i_ d df 

dz; dz /•] 

+ T 1 -' 0 ~d( +IJ w\ k ,5b) 

The second derivatives are 

d?v did du d?(i . , , 

r “ = [7J y + ' 2 ^ {y f- Ui) + " 77 {y f - pi) \ J 

d?v „ n did du, . . d?n , * / n \ 

+ [77 z 7 + u 77 7 ' j k (6a) 
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With equations (4) through (6) all the geometric 
parameters required for the surface wave calculations 
can be computed. 

4.3. Ray Tubes and Ray Strips 

The surface wave speed is assumed to be constant. 
Consequently the ray paths, which by Fermat's prin- 
ciple are least time paths, will also be paths of min- 
imum distance, or geodesic lines. Given a point on 
the surface and an initial direction, a geodesic hav- 
ing these initial conditions can be calculated by a 
marching procedure. (See ref. 4 for details.) 

For the purpose of analyzing the surface rays, 
two surface coordinate systems are utilized. One is 
the natural (u,v) coordinate system that is used to 
define the surface itself. For this system, the metric 
coefficients are 

E=r u -r u (7a) 

F=r u -r v (7b) 

G=r v -r v (7c) 

and the discriminant is 



Figure 3. basic geometry illustrating surface parameters 
u, v, t, and i). 

da = \fEx- = c s dt (11) 

The incremental arc length normal to the rays (that 
is, along the surface wave fronts) is 

dn = \J G r dr) ( 12) 

where G r is to be determined and dr) is the arbitrary 
inclement that distinguishes two adjacent rays. It 
can be taken as the distance between the rays along 
the initial excitation line. The metric coefficient 
G r is determined as follows. As each point (u r ,v r ) 
is computed, the local values of e, /, g, and H 
are computed. From these, the local value of the 
Gaussian curvature K is calculated by the formula 


H = \/EG - F' 2 (8) 


K = 


~ f 2 

H 2 


(13) 


The local surface normal is 

N = (9) 

and the curvature coefficients, or coefficients of the 
second fundamental form, are 


e — r uu 

■ N 

(10a) 

f = I’m: 

■ N 

(10b) 

g = r vv 

■N 

(10c) 


where the notation is consistent with that of refer- 
ence 9. 

The second coordinate system used consists of 
time t along a family of rays all initiated at a given 
time t 0 and a surface coordinate normal to the rays 
(fig. 3). Thus, incremental distance along the rays is 


Then G r is computed by numerical integration of the 
equation 

^E + c !av 57= o (i4) 

(The factor c 2 was inadvertently omitted in equa- 
tions (25) and (29) of ref. 4.) The geodesic curvature 
of the surface wave front is also needed. It is obtained 
from G r by the formula 


K g,n — 


1 d\J Gr 


d\/U7 


\J G r da fy \/ c, r dt 


(15) 


The geodesic curvature of the ray itself is zero by 
definition. The local surface normal curvature n 
in the direction of the ray is 


n,r 


ihere 


_ e + 2/A + g\ 2 
hn ' r E+2FA + GX 2 


dv 

du 


(16) 

(17) 
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on the ray. The normal curvature in the direction of 
the wave front is similarly 


_ c + 2/A„ + g 
Kn ' n ~ e + 2F\ n + ( /a;; 


(18) 


where X n is the direction orthogonal to the rays. 

The initial locations and directions for rays ex- 
cited by an off-surface source are determined as fol- 
lows. Let I = I\i + I- 2 J + Isk denote the direction 
vector from the source to an arbitrary point (u, v) on 
the surface. Then at (u, v), the conditions 


I(u,v) ■ N(u,v ) = cos 9 (19) 

where 

sin0= — (20) 

c s 

must be satisfied. The locus of the initial excita- 
tion line is then found by numerical solution of equa- 
tion (19). In practice, one usually selects values of 
u or v and solves equation (19) for the remaining 
variable (e.g., search along a u = Constant line for a 
value of v that satisfies the equation). For a far-held 
source, I is constant, and for Franz (creeping) waves 
initiated by grazing rays, cos 0 = 0. 

To determine the initial ray direction A 0 , we take 
a step dr from r 0 at (u, v): 


Let A s denote the amplitude of the surface wave. 
Then, following the analysis of reference 2 (eqs. (4) 
through (6)), we assume that the energy radiated 
into the medium from an increment da of a ray strip 
is proportional to A? s and to the area da dn of the 
ray strip that radiates this energy. Thus, 


Ps c s 


(<7 + da ) dn (a + da) — flj( a) dn da 


— Ps c s 


2p,(a)A^ dn da 


(23) 


where, again consistent with reference 2, 2p.(a) is 
taken to be the proportionality constant. The at- 
tenuation factor p. is preferably determined by ex- 
periment if such data are available, but otherwise, 
it is determined by solution of the wave equation 
for a canonical problem. Expressions for p. are tab- 
ulated in reference 2 for several types of boundary 
conditions. 

Equation (23) leads to the differential equation 

^(/4g dn p s c s ) = -2 p(A 2 s dn p s c s ) (24a) 
d(A 2 dn p s c s ) = —2 p(A^ dn p s c s ) da (24b) 
whose solution is 


dr = (x u du + x v dv)i + (y u du + y v dv)j 


+ (z u du + z v dv)k 


( 21 ) 


j 4.s(cr) — A up 



1/2 

exp 



p(a) da 


(25) 


If the surface rays are launched by grazing rays 
(0 = 90°), dr must be parallel to I. If 9 < 90°, dr 
must still be parallel to the projection of I onto the 
tangent plane at r 0 . In either case, the components of 
dr must be proportional to those of I|| = I— (I ■ N)N . 
For example, 


yu du + y v d» _ Ip 
x u du + x v dv I|| 1 

V u + _ % 

X U “1“ X V *^0 I|| ? l 


(22a) 

(22b) 


A 


o 



x v y v 


(22c) 


Clearly, any two components of dr could have been 
chosen for this calculation. 


The phase on the ray strip is 

if) = <f) 0 + ka ( 26a) 

<t> = <f> - V kc s t (26b) 

For the rays radiated from the surface, the GTD 
theory assumes that the amplitude of the radiated 
ray is proportional to the amplitude of the surface 
ray (ref. 3): 

A r = RA 8 (27) 

For the radiated ray tube, the energy flux is 

A.E r = A r dr p m c m (28) 

where dr is the cross-sectional area of the ray tube. 
At the surface (fig. 4fe 
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(33b) 



Figure 4. Geometry of ray tube cross-sectional ai'ea. 


dr s = dn da cos 


(29) 

Equating the expression in equation (28) to the mag- 
nitude of the energy flux in equation (24b) and using 
Equations (27) and (29) yield 

2pA 2 p s c s da dn = R 2 A 2 p m c m dn da cos 9 (30) 

from which 


R 2 = 


Ps^S 

cos 9 


(31) 


This method of obtaining the radiation coeffic ient 
from the attenuation factor p could also be applied 
inversely. That is, if R were determined directly from 
a solution to a canonical problem, then p could be 
obtained from R. 

Let p denote the direction vector of a radiated ray 

p = sin 9 T + cos 9 N ( 32 ) 

and let q denote the unit vector normal to p in the 
T, N plane. (See fig. 5.) The change in p over the 



Figure 5. Orientation of unit vectors in tangent (Tpu.) plane 
and in radiation (T, N) plane. 


q ■ ) da = ( cos 9 T — sin 9 N ) ■ (n r sin 9 N 

da 


— n r cos 9 T ] da 


(33a) 


q ■ \ da = k t da 

da J 


Let K n denote the curvature vector of the surface 
wave front. The variation of p over the distance dn 
along the surface wave front is 

(Kn ■ p) = (kii.h, cos 9 + n n g sin 9) dn = K p dn 

(34) 

The angular increment along the ray over distance 
da is given by equation (34), and the angular incre- 
ment over dn normal to the ray is given by equa- 
tion (33b). Thus the spreading of rays emitted from 
this incremental area is sn r da in the N, T plane, and 
SKp dn normal to this plane. The actual ray tube 
area at distance s from the surface, measured along 
the ray, is obtained by adding these increments to 
the increments at the surface (eq. (29)), which yields 
(see fig. 4) 


dr s = (cos 9 + sK r ) da (l + k p s) dn 


(35) 


By formulas (27), (29), and (35), the amplitude 
A p at point P is given by 


A' 2 (cos 9 da -\- K r s da)(dn + n p s dn) 


A% 


da dn cos 9 


(36a) 


A 2 


— = ( cos 9 + k 7 ,s)(1 + K p s) sec 9 (36b) 


A 


where s is the distance from the radiation point P r 
to P. A more general form of equation (35) is given 
in reference 9. 

A problem arises with Franz waves for which 
9 = f-, since in this case R in equation (31) increases 
indefinitely. The singular nature of the radiation 
coefficient and the caustic condition at the surface 
are discussed at length in the literature (refs. 2 and 3, 
for example). However, this singular character has 
a somewhat artificial quality, since it arises from 
the manner in which the physical phenomena are 
modeled rather than the actual physical process. If 
we relate the amplitude at P directly to the surface 
wave amplitude by equations (27), (31), and (36), we 
get the nonsingular result 


4 2 

—A 

A% 


A P \ 2 ( A . r ' 2 


.4 




= 


cos 6 


(cos 9 .s /> ]( I p /.' pS p ) 

2 P feG 

(cos 6 + K r Sp)(] + Kp .Sp) p m c,- 


(37) 
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If the surface wave is excited by an off-surface 
source at Q (fig. 6), the initial amplitude on the 
surface at the excitation point Q; is related to that 
of the incoming wave by an excitation coefficient e: 



A s ,o — e(Qi) A ; i (38) 


path. Furthermore, the attenuation due to spread- 
ing is clearly the same for both the forward and the 
reverse directions. Consequently, when the forward 
and reverse amplitudes obtained from equation (40b) 
are equated, the reciprocity principle yields 



Sq (cos 9 + K r S )( 1 + Kp-s) 


Sp [(cos 9 + K r S )( 1 + KpS)] Q 

Thus, in order for reciprocity to be satisfied (ep = 
£q), we must have 


£ Q 


f Q 

(cos 9 + K r s)( 1 + KpS) Q 


(42) 


In this formula, (np)Q. is the component in the 
direction —I of the curvature vector of the initial 
wave front, which is found by solving equation (19) 
numerically, as was discussed earlier. Also, K r is 
the surface curvature in the surface direction A 0 as 
computed from equations (16) and (22c). 

The phase u it along the ray path is 


where A- is related to a source strength A at Q by 


,2 A 2 

M = — 

S Q 


(39) 


(See fig. 6.) Combining equations (38), (37), (25), 
and (36) yields 


A 
' Ip 



(40a) 


2 // 

(COS 0 “t“ AV-Sp)(l “t“ KpS i . j 


J P r 



x exp 



(40b) 


where the subscript P r denotes that the quantity is 
evaluated at the radiation point P r . 

To determine the excitation coefficient we ap- 
ply the reciprocity principle, which states that an 
equal field is obtained at Q if the source is located at 
P instead of at Q. Since for both cases the surface ray 
path is the same, the spreading factor dn Q /dn is the 
same if the directions are reversed on the path. This 
follows from equations (12) and (14), which demon- 
strate that the ray strip spreading is determined by 
the distribution of the Gaussian curvature along the 


(j) = LO 



+ ^+-[ F 

c m c s J Qj 



(43) 


Appropriate phase jumps of integer multiples of 7t/2 
must be included if the ray traverses caustic points. 

4.4. D iscussion of Related Problems 

In deriving the relationship between the attenua- 
tion factor and the radiation coefficient, the attenu- 
ation is assumed to be caused only by the radiation 
of energy from the surfaces If the waves penetrate 
the surface and attenuation also occurs because of 
dissipation or inhonrogeneity, such attenuation must 
be included in the analysis. Reference 9 considers 
waves in a ray tube adjacent to the surface and ob- 
tains a somewhat different result from that of Levy 
and Keller (ref. 2) for waves satisfying the impedance 
boundary condition at the surface. For certain elas- 
tic waves, a change in wave type occurs as the wave 
in the medium excites a surface wave or is emitted 
from the surface (ref. 10). As might be expected, in 
this case the factor p m c m / p s c s in equation (31) is not 
unity. 

As was mentioned earlier, vector waves are often 
treated by synthesizing the solutions from those of 
scalar waves. The polarization of the waves, however, 
depends on the torsion of the ray path. The torsion 
may be computed at each point on the path by first 
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computing N and the quotient dN /da and using the 
formula (ref. 11, p. 201) 


T(<r) = T xN ■ — (44) 

da 

Actually, the quantity dN / da is the crucial factor 
since, for perfectly conducting surfaces, the wave is 
polarized normal to the surface. 

4.5. Wave Fronts and Strengths 

In section 4.3 it was shown that one can compute 
surface geodesic ray paths and express the location of 
points on these rays in terms of coordinates t and i) f , 
where t is time elapsed from the point of excitation 
and i] distinguishes the individual rays. Now the 
local wave front surface at an arbitrary time t* is 
expressed in the vector form of equation (1) in order 
to more easily study the surface geometry of the 
wave front. For this purpose the surface variables 
u, v in equation (1) are taken to be ( t,i ]). Since 
c s is assumed to be constant, we could replace the 
coordinate t with a = c. s U Denoting a point on a 
surface ray by x(f,?/), the ray emitted from x(f,?/) 
passes at time t* through the point 

y (t , rj) = *(t, rj) + c m(t* ~ t) p(t, i)) (45) 


where, from equation (32) and the Frenet formulas, 


dp 

da 


d.T dN 

— + cos 9 — — 

da da 


(49 a) 


dp 

da 

dp 

da 


sin 9 K r N + cos 
K r q cos 9T r u. 


^K r T + T r uj (49b) 
(49 c) 


w 1 1 qt i T,. is the torsion of the surface ray and N 
points are in the opposite direction from the curvi- 
linear vector of the geodesic ray path. 

Substituting equations (47) and (49c) into (48b) 
yields 


— = e^sin 9 p+ c s cos 9 q — c m p + sc^Kj.q 

ill 



+ sc s cos 9 T j.u 


(50a) 

which, 

becomes 

with the 

assumption 

of equation 

to 

o 

II 

: Cg ( COS 9 

+ SKj.)q + Cg s 

cos 9 T j.u 

(50b) 


where p is the direction vector defined by equa- 
tion (32). Utilizing equation (45), which has the 
required form of equation (1), we can derive an ex- 
pression for the wave front surface area element and 
thereby verify formula (35) for the ray tube cross- 
sectional area. 


plane by q (fig. 5). Then 

q(t,rj) = cos 9 T — sin 9 N 
We can solve equations (32) and (46) for T 
T = sin 9 p + cos 9 q 


p,q,u form an orthonormal set. file 
y (t,rj) with respect to t is 

dy dx dp 

W = df~ CmP + s li 

This equation can also be written 

d T _ ^ . dp 

i, — c .sd c mP + sc .s i 
dt da 


Differentiating y with respect to i] yields 

dy dx dp 

dr) dr/ dr) 

where 


dp 

dn 

( . dT 

sm 9 — 

V dn 



dp 


dr) 

(46) 

dp 


dr) 

(47) 

dp 


dr) 

■ N,T 

T hen 

dp 

;ive of 

dr) 


Thus 

(48a) 



The wave 

(48b) 



?-=Vg; k 


= Vg7( i + sk p )i 


(51) 


(52a) 


dN \ 

W 

? Kn,n'u) (52c) 
(5 2d) 

(53) 


dt dr) , 


(54) 
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which yields (with eqs. (11), (12), (50b), and (53)) 


P Q 


u 


AA = 


(1 c s (c os (9 + k,,.s) c s cos Os T,, 

0 0 '/( o* ' 1 + SK.p) 


da dn 

c s \/G~ 


(55a) 


= (l + SK^)(cos 0 + K r s) dn dap (55b) 

Thus p is the wave front normal direction, and the 
ray tube cross-section area is 

dr = (1 T s cos 9 + n r s) dn da (56) 

which is consistent with that determined in the pre- 
vious section (see eq. (35)) and, thus, provides an 
independent verification of the theory. One could 
also differentiate equations (50b) and (53) to obtain 
the curvatures of the wave front surface in terms of 
those of the radiating surface. 

It is of interest to study the wave front equation 
(eq. (45)) in the limiting two-dimensional case. That 
is, we consider waves radiated from a cylinder hav- 
ing a smooth cross-section shape which is arbitrary 
within the constraint that none of the radiated rays 
are incident on the surface. Then equation (45) be- 
comes 

y (t) = x(f) + c m (t* - t)p(t) (57) 

and equation (50b) simplifies to 

d f = (C, CCS 6 + K r )q(t) (58) 


In particular, if we treat Franz (creeping) waves 
initiated by a ray grazing a convex cylinder at a 
point Q (fig. 7), then, since for this case c s = c m , 
the constant phase at the wave front corresponds to 
a constant distance the rays travel over the surface 
from Q to x(f) and from x(f) to y. Consequently, 



Figure 7. Shape of diffracted wave front emitted by Franz 
waves. 


the diffracted wave front represents the involute of 
the curve that defines the cross-section shape of the 
diffracting surface. In this case, equations (57) and 
(58) become, respectively, 

y (t) = xi/i + c m (t* - /!/ (/! (59) 

and 

% = Kr(t)N(t) (60) 


5. Sample Calculations 

In the following examples, the ray calculations 
were performed by the method of reference 4. The 
high accuracy of the method was verified by compar- 
ing a numerical calculation with the known helical 
geodesic on a circular cylinder. For the illustrative 
example shown in figure 8, the surface was defined 
with the initial cross section taken to be an ellipse 
having a ratio of minor to major axis of 1/3, and the 
base section is an ellipse having a ratio of 0.58. This 
surface is smooth and convex but not separable. Fig- 
ure 8(a) shows the results of a Franz wave calculation 
with an illuminating beam having direction numbers 
(1,7,2). For clarity in presentation, only a portion 
of the surface wave and the corresponding diffracted 
wave are shown. Figures 8(b) and (c) show a simi- 
lar calculation from different viewing angles but with 
the surface wave speed c s three times greater than the 
wave speed c m in the medium. 


6. Resonant Phase Matching 

If we compute a geodesic ray path that completely 
circumnavigates a smooth body, the path may close 
on itself. That is, it may return to the initial point 
propagating in the initial direction. This always 
occurs on spheres and in two-dimensional situations. 
Furthermore, it occurs for propagation along the 
maximum or minimum meridian line on a spheroid 
or any smooth closed surface that has a maximum or 
minimum meridian. 

When such a closed path exists, a resonant condi- 
tion occurs for a frequency such that the phase accu- 
mulation over the closed path is an integer multiple 
m of 27T . Since surface waves over closed surfaces 
form caustics, this phase accumulation must include 
a shift of 7r/2for each time the ray traverses a caustic 
point. In equation form 


Ul 

Cs 


ds — m — = 2 7T777. 

2 


(61) 
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Surface rays 


(c) 9 = cos 1 3; view from above surface. 
Figure 8. Concluded. 



where m is the number of caustic traversals. The 
calculation of resonant condition rays for separable 
surfaces and the relation of the corresponding res- 
onant frequencies of the natural frequencies of the 
body have been reported by Uberall and his asso- 
ciates in several papers (refs. 12 and 13). In order 
to determine a resonant condition for a nonsepara- 
ble surface, one must first find a closed surface ray 
path and then determine the number of times the ray 
traverses a caustic over this path. 

The body used for the surface calculations of fig- 
ure 8 was investigated to locate closed ray paths. 
Obvious closed paths exist along the maximum and 
minimum meridian lines. In addition to these cases, 
another type of closed path was found (by an itera- 
tive procedure), as shown in figure 9. Here, the ray 
circumnavigates the body twice before the closed cir- 
cuit is complete. The caustic points on this path are 
points where G r vanishes. For illustrative purposes, 


an adjacent ray was calculated in order to display the 
ray strip with its two caustics in the closed ray path. 

Thus, for this case, if a c is the closed path length, 
resonance occurs for 

Wrn — —(2m + 1)7T (62) 

In this situation, if the wave is a Franz wave, the 
effect of the resonance is seen in the rays diffracted 
from the surface along this closed ray path. How- 
ever, if the surface wave is initiated by a grazing 
ray from a far-held source, the effect is not observed 
at the source because no ray is diffracted along this 
path in the direction of the source. It could be de- 
tected by a detector (microphone or antenna) in a 
bistatic arrangement. This kind of closed ray path 
might also affect the mutual interference of surface 
mounted antennas, although the primary concern for 
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this problem is the direct geodesic distance between 
the antennas. 

As was mentioned, several papers by U be rail and 
his associates (refs. 12 and 13) have discussed the 
relationship between surface wave frequencies and 
the normal modes of an elastic body, and calculations 
were performed for some separable surfaces. With 
the present capability to determine ray paths on a 
broad class of nonseparable surfaces, the potential 
exists for a variety of applications relating to elastic 
phenomena. Consider, for example, Lamb waves in 
a thin shell structure (with the usual restriction that 
the radii of curvature are large relative to the wave 
length, so that the wave speed is essentially constant 
for a given frequency). Then, in addition to the usual 
modes associated with ray paths along the maximum 
and minimum meridian lines, the possibility exists for 
a weaker kind of resonance associated with closed ray 
paths such as that shown in figure 9. 

7. Concluding Remarks 

A theory has been presented for computing sur- 
face waves and waves radiated from a smooth surface. 
With the assumption that attenuation of the surface 
wave is due only to radiation, the radiation coefficient 
was derived in terms of the attenuation factor. Then 
the excitation coefficient was determined by the reci- 
procity condition. Formulas for the shapes and the 
spreading of the radiated waves were derived, and 
some sample calculations were given. An investiga- 
tion of resonant phase matching for nonseparable sur- 
faces was presented with a sample calculation. Also 
included was a discussion of how such calculations 
might be related to resonant frequencies of nonsepa- 
rable thin shell structures. A description was given 
of nonseparable surfaces that can be modeled in the 
vector form that facilitates use of the appropriate 
formulas of differential geometry. 

NASA Langley Research Center 
Hampton , VA* 23665-5225 
February 20, 1092 
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